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ABSTRACT 

Bulges are commonly believed to form in the dynamical violence of galaxy collisions and mergers. Here 
we model the stellar kinematics of the Bulge Radial Velocity Assay {BRAVA), and find no sign that the Milky 
Way contains a classical bulge formed by scrambling pre-existing disks of stars in major mergers. Rather, the 
bulge appears to be a bar, seen somewhat end-on, as hinted from its asymmetric boxy shape. We construct a 
simple but realistic A^-body model of the Galaxy that self-consistently develops a bar The bar immediately 
buckles and thickens in the vertical direction. As seen from the Sun, the result resembles the boxy bulge of 
our Galaxy. The model fits the BRAVA stellar kinematic data covering the whole bulge strikingly well with no 
need for a merger-made classical bulge. The bar in our best fit model has a half-length of ^ 4 kpc and extends 
20° from the Sun-Galactic Center line. We use the new kinematic constraints to show that any classical bulge 
contribution cannot be larger than ^ 8% of the disk mass. Thus the Galactic bulge is a part of the disk and not 
a separate component made in a prior merger Giant, pure-disk galaxies like our own present a major challenge 
to the standard picture in which galaxy formation is dominated by hierarchical clustering and galaxy mergers. 

Subject headings: Galaxy: bulge - Galaxy: kinematics and dynamics - galaxies: kinematics and dynamics 



1. INTRODUCTION 

Astronomers commonly make the Copernican assumption 
that our Milky Way is in no way unusual; then they can 
exploit the fact that we live in it to study galaxy forma- 
tion in special detail. Past assumptions about our Galactic 
bulge grew out of our developing understanding of galaxy 
formation. It is well known that spiral galaxies consist of 
three main components, an invisible dark matter halo, an 
embedded, flat disk, and a central bulge. The bulge of our 
Galaxy is >99% ma de of stars that are at least 5 Gyr old 
( Clarkson et alj|2008l) with a wide range of metal abundances 
JMcWiUiam & R ich 1994; Ful bright et al .^20Q6: Zocca U etal] 
120081) . In this respect and many others, big bulges are sim- 
ilar to (diskless) elliptical galaxies. The formation of el- 
lipticals is well understood. Hierarchical gravitational clus- 
tering of initial fluctuations in the cosmological density re- 
sults in galaxy collisions and mergers that scramble flat disks 
into rounder elliptic als (Toomre 1977; White & Rees 197^ 
ISteinmetz & Navarro 2002; Nakasato & Nomoto 2003). Sig- 
nificant energy has been invested in developing this very suc- 
cessful theory of galaxy formation, and it was natural to think 
that our Galactic bulge is a product of it. There is little danger 
that the picture is fundamentally wrong (Binnev 2004). 

But it is incomplete. The theme of this paper is that our 
Galactic bulge is indeed normal but that it is prototypical of 
different formation processes than are usually assumed. A 
complementary suite of evolution processes shapes isolated 
galaxies. They evolve by rearranging energy and angular mo- 
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mentum; this grows central components that masquerade as 
classical bulges but that formed directly out of disks with- 
out a ny collisions (Kormendy 1993; Kormen dy & Kennicutj 
120041) . To distinguish them from merger remnants, we call 
them "pseudobulges". They come in two varieties. Some 
are flattened; they are grown out of disk gas transported 
inward by nonaxisymmetries such as bars. Another va- 
riety is recognized only in edge-on galaxies. When bars 
form out of disks, they buckle vertically and heat them- 
selves into thi ckened structures that look box-shaped when 
seen edge-on (ICombes & Sandersl[l981l:ICombes et al.lll99"ol: 
iRaha et anil991l) . Our Galaxy contains suc h a box-shaped 
(lMaiharaetai]I j978; Weiland et al. 1994; D weket alJ[T995h 
pseudobulge (Korme ndv & Kennicuttil200 4). 

The identification of this boxy structure as an edge- 
on bar is particularly compelling beca use infrared irnagery 
shows a parallelogra r n-shaped distortion (iMaihara et al.l 19781 : 
IWeilandetaIlll994l: iDweket alj|1995h fliat is naturally ex- 
plained as a perspective effect: the near end of the bar is closer 
to us than the far end. So its vertic al extent is taller on the 
near s ide than on the far side (Blit z & Spergell[l991h . IZhaol 
(Il996h developed the first rapidly rotating bar model that fitted 
this di stortion. Zhao's model was based on the lSchwarzschildl 
(119791) orbit superposition technique, so it was self-consistent 
and in steady state, but it did not evolve into that state from 
plausible initial conditions. Also, little stellar kinematic data 
were available to constrain Zhao's ste ady-state mod el and 
eai-ly A^-body models (Fux 1997, 1999; S evenster et al.lfl9 99). 
and subsequent radial velocity data from a survey of plane- 
tary neb ulae, although com pared with a range of dynamical 
models dBeaulieu et alj|200 0). led to only limited conclusions 
because of the small numbers and uncertain population mem- 
bership of the planetary nebulae. 

In this paper, we simulate numerically the self-consistent 
formation of a bar that buckles naturally into a thickened 
state, and we scale that model to fit new kinematic data on 
bulge rotation and random velocities. The radial velocity 
observations are provided by the Bulge Radial Velocity As- 
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Fig. 1. — Upper three panels: Face-on and side-on views of the surface 
density of our best-fitting model as seen from far away. The Sun's position 
8.5 kpc from the Galactic center is marked along the +.y axis. The Galaxy 
rotates clockwise as seen in the face-on projection. Bottom panel: Model 
surface brightness map in Galactic coordinates as seen from the Sun's loca- 
tion. Our perspective makes the box-shaped, edge-on bar look taller on its 
nearer side. The Galactic boxy bulge is observed to be similarly distorted. 



say {BRAVA; iRich et alj|2007t [Howard et al]|2008h . This is 
a spectroscopic survey of the stellar radial velocities of M- 
type giant stars whose population membership in the bulge 
is well established. These giants provide most of the 2 /^m 
radiation whose box-shaped light distribution motivates bar 
models. BRAVA emphasizes measurements in two strips at 
latitude b = -4° and b = -8° and at longitude -10° < / < 
-1-10°. A strip along the minor axis (/ = 0°) has also been 
observed. We use nearly 5,000 stellar radial velocities in this 
report. A pr eliminary analysis o f data found strong cylindri- 
cal rotation ([Howard et al.ll2009h consistent with an edge-on, 
bar-like pseudobulge, although a precise fit of a bar model to 
the data was not available. This success leads us here to con- 
struct a full evolutionary A^-body model that we can fit to the 
radial velocity data. 

2. MODEL SETUP 

We use a cylin dr ical particle-m esh ^code 

(iSellwood & Valluril UqQTI: iShen & SellwoodI HOO?) to 
build fully self-consistent A'-body galaxies. It is well suited 
to study the evolution of disk galaxies: we model the disk 
with at least 1 million particles to provide high particle 
resolution near the center where the density is high. We 
try to construct the simplest self-consistent A^-body models 
that fit the BRAVA data, avoiding contrived models with 



too many free parameters. Initially, they contained only an 
unbarred disk and a dark halo. The profile of the Galactic 
halo is poorly constrained observationally; we adopt a rigid 

pseudo-isothermal halo potential $ = jV^ ln(l -I- ^). Here Vc 

^ 250 km s"' is the asymptotic circular-orbit rotation velocity 
at infinity, and R^ = 15 kpc is the core radius inside which 
the potential is effectively constant. This halo gives a nearly 
flat rotation curve of ~ 220 km s"' between 5 to 20 kpc. A 
simple halo form allows us to run many simulations quickly; 
this is important for a parameter search such as the present 
one. A rigid halo also omits dynamical friction on the bar, but 
the central density of the cored halo we adop t is low enough 
so that friction will be very mild (e.g. Debatt ista & SellwoodI 
I2000t lAthanassoula & MisiriotisI 12002). More importantly, 
we are mainly interested in the bulge, which is embedded 
well interior to R^. So the exact profile of the dark halo at 
large radii is not critical. We will explore more sophisticated 
halos in a future study. 
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Fig. 2. — (top): Mean velocity and velocity dispersion profiles of the best- 
fitting model (black lines) compared to all available kinematic observations. 
The left two panels are for the Galactic latitude h = —4° strip; the middle two 
panels are for the b = -8° strip; and the right two panels are for the / = 0° 
minor axis. The black diamonds and their error bars are the BRAVA data; the 
green diamonds are forM-type giant stars (Rangwala et al. 2009), and the red 
triangles are the data on red clump giant stars (Rangwala et al. 2009). This is 
the first time that a single dynamical model has been compared with data of 
such quality. The agreement is striking. 



3. RESULTS AND DISCUSSIONS 

In our models, a bar develops self-consistently from the 
initially unbarred, thin disk. Bar formation enhances the ra- 
dial streaming motions of disk particles, so the radial velocity 
dispersion quickly grows much bigger than the vertical one. 
Consequently the disk buckles vertically out of the plane like a 
fire hose; this is the well known buckling or corrugation insta- 
bihty (Toomre 1966; Combes et al. 1990; Raha et al. 1991). It 
raises the vertical velocity dispersion and increases the bar's 
thickness. This happens on a short dynamical timescale and 
saturates in a few hundred million years. The central part of 
the buckled bar is elevated well above the disk mid-plane and 
resembles the peanut morphology of many bulges includin g 
the one in our Galaxy ("Combe s etaPl 19901: iRaha et al.lll991h . 

Out of a large set of A^-body models, we find the one that 
best matches our BRAVA kinematic data after suitable mass 
scaling. The barred disk evolved from a thin exponential disk 
that contains Md = 4.25 x lO'^M©, about 55 % of the total 
mass at the truncation radius (5 scale-lengths). The scale- 
length and scale-height of the initial disk are ~ 1 .9 kpc and 
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Fig. 3. — Model results compared to the data from the b = —4° major-axis 
strip as we vary the bar angle relative to the line that connects the Sun to the 
Galactic center. The red, black, green, blue, and cyan curves are the model 
results for bar angles of 0°, 20° (the adopted best-fit value), 45°, 67.5°, and 
90° , respectively. All data points are plotted as black diamonds regardless of 
source. Clearly bar angles of 0° and 20° are favored over larger values. 

0.2 kpc, respectively. The disk is rotationally supported and 
has a Toomre-Q of 1.2. The amplitude of the final bar is in- 
termediate between the weakest and strongest bars observed 
in galaxies. The bar's minor-to-major axial ratio is about 0.5 
to 0.6, and its half-length is ^ 4 kpc. Figure 1 (top three pan- 
els) shows face-on and side-on views of the projected density 
of the best-fitting model. A distinctly peanut shaped bulge is 
apparent in the edge-on projection. Figure 1 (bottom panel) 
shows the surface brightness distribution in Galactic coordi- 
nates as seen from the Sun's vantage point. Nearby disk stars 
dilute the peanut shape, but the bar still looks boxy. More- 
over, from close up, an asymmetry in the longitudinal direc- 
tion is apparent; this means that the bar cannot be aligned with 
the direction from the Sun to the Galactic center Rather, its 
near end is at positive Galactic longitude, so it looks taller in 
that quadrant, and it extends farther from the Galactic cen- 
ter on the near side than on the far side. Both the boxy 
shape and the asymmetry are in good agreement with the mor- 
phology revealed by the COBE satellite near-infrared images 
dWeiland et al..l994; Dweket al..l9 95). 

Figure 2 compares the best-fitting model kinematics (solid 
lines) with the mean velocity and v elocity disp ersion data 
from the BRAVA and other surveys (iRangwala e t al. 2009). 
All velocities presented here have been converted to Galac- 
tocentric values (the line-of-sight velocity that would be ob- 
served by a stationary observer at the Sun's position). For 
the first time, our model is able simultaneously to match the 
mean velocities and velocity dispersions along two Galactic 
latitudes (-4° and -8°) and along the minor axis. 

Figure 3 constrains the angle between the bar and the line 



that connects the Sun to the Galactic center. It compares the 
model results with the data in the b = -4° major-axis strip 
as we vary the above angle. Clearly the smallest bar angles 
give the best match to the velocity dispersions. Intriguingly, 
we find that the velocity dispersions provide much stronger 
constraints than the mean velocity profile. A bar angle of 
0° also matches the kinematics well. However, the photo- 
metric asymmetry excludes a bar that is pointed at the Sun. 
We therefore conclude that the overall best-fitting model has 
a bar angle of ^ 20°. Other studies converged on a similar 
bar angle (Staneket al. 1997; Fr eudenreich 1998; Fux 199'^ 
ll999tlBissantz & Gerhardll2002l) . The excellent match to the 
data in Figures 1-3 strongly supports the suggestion that 
the boxy pseudobulge of the Milky Way is an edge-on, buck- 
led bar that evolved from a cold, massive disk. The thick- 
ened disk in the pseudobulge-forming process may have con- 
tributed to the thick disk of the Milk Way, as hinted from 
chem ical similarities of Galactic bulge and loca l thick disk 
stars jAlves-Brito et al.ll2010l;lBensbv et al.ll2010l) . 

The model in Figures 1-3 contains no classical bulge 
component. Could a small classical bulge also be present? 
Could it have been spun up by the formation of a bar, flat- 
tened thereby and made hard to detect? To constrain such 
multi-component models with BRAVA kinematics, we also 
constructed models with a pre-existing classical bulge. The 
distribution function for the hve classical bulge c omponent 
was generated iteratively ( Debattista & Sellwoodl200Qi) to en- 
sure that both the disk and the classical bulge were initially in 
equilibrium, and we required that the bulge parameters are 
clos e to the fundamental pla ne for classical bulges and ellipti- 
cals (iKormendv et al.ll2009l) . The setup of the disk is the same 
as in the disk-only model. We show in Figure 4 that inclusion 
of a classical bulge - one widely thought to be typical of Sbc 
spiral galaxies like our own - greatly worsens the model fit 
to the data. The degradation is especially obvious along the 
Galaxy's minor axis. Including a classical bulge with just 8% 
of the disk mass considerably worsens the fit of the model to 
the data. Our models rule out that the Milky Way has a sig- 
nificant classical bulge whose mass is >~' 15 % of the disk 
mass. 

Could a smaller, merger-built bulge hide inside the boxy 
bar? The only result that we are aware of that might point 
to such a conclusion is the observed drop in stellar metal 
abundances with increasing h eight above the Galactic p lane 
JZoccaU et a l. 2008; Zoccali 2()l3). IZoccali et alj (l2008l) ar- 
gue that this means that the bulge must consist of both a clas- 
sical and an edge-on bar component. However, no kinematic 
gradient or transition corresponding to the abundance gradi- 
ent is observed. Moreover, an abundance gradient can be 
produced within the context of secular pseudobulge forma- 
tion if some of the vertical thickening is produced by resonant 
heating of stars that scatter off the bar (Pfenniger & Normao 
Il990h . If the most metal-poor stars are also the oldest stars, 
then they have been scattered for the longest time and now 
reach the greatest heights. 

Our results have important implications for galaxy forma- 
tion. We demonstrate that the boxy pseudobulge is not a 
separate component of the Galaxy but rather is an edge-on 
bar. Bars are parts of disks. To be sure, the stars in our 
Galactic bar are older than most disk stars. But those stars 
could have formed over a shor t period of t i me but lo ng be- 
fore the bar str ucture forme d dWysd [19991; iFreemanI [200 8), 
their old age (Zocc ali et all 120031; iFulbright et all 120071) is 
therefore not an argument against the internal secular evolu- 
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Fig. 4. — Fits to the kinematic data (cf. Figure 2) of models that include a pre-existing classical bulge. The heavy black lines from Figure 2 represent the model 
without a classical bulge. The red, green, and blue lines are for models whose classical bulges have masses of 8%, 15%, and 30%, respectively, of the disk mass 
Miisk- Including a classical bulge significantly worsens the model fits to the data, especially along the minor axis. 



tion model. Our kinematic observations show no sign that 
the Galaxy contains a significant merger-made, "classical" 
bulge. So, from a galaxy formation point of view, we live 
in a pure-disk galaxy. Our Galaxy is not unusual: it is very 
similar to another giant edge- on galaxy with a boxy bulge, 
NGC 4565. lKormendy & B^ ntine (2010) recently show that 
NGC 4565 does not contain even a small classical bulge com- 
ponent and that it therefore is another giant, pure-disk galaxy 
that contains no sign of a merger remnant. In fact, giant, pure- 
disk galaxies are common in environ ments like our own tha t 
are far from rich clusters of galaxies (iKormendv et al.ll2010l) . 
Classical-bulge-less, pure-disk galaxies present an acute chal- 
lenge to the current picture of galaxy formation in a Universe 
dominated by cold dark matter - growing a giant galaxy via 



hierarchical clustering (Vc ~ 220 km s"' in the Milky Way) 
involves so many mergers that it seems almost impossible to 
avoid forming a substant ial classical bulge (Peebles & Nusse^ 
l2010l:lAgertz et alll2010l) . How did our Galaxy grow so large 
with no observational sign that it suffered a major merger af- 
ter the time 9-10 Gyr ago (e.g. Winget & Kepler 2008) when 
the first disk stars formed? 
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